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The variational method and linear combination of atomic orbitals (LCAOs) 

Suppose we wanted to solve the eigenvalue problem Hip = Eip for an n-atom system, but we have no idea 
what the mathematical expression for wavefunction ip looks like. The variational method is a tool that 
nevertheless allows us to gain insight into its properties. First, recall that if ip is a normalized solution to 
the eigenvalue problem, the expectation value of energy is given by 


E = E 


ip* ip 


ip* (Eip) 


ip* Hip. 


(All integrals in this handout are over all space ]R 3 , with respect to volume element dr = dxdydz.) The 
variational theorem allows us to make and optimize a guess. 

Let Eq be the smallest eigenvalue (i.e., ground state energy) with corresponding wavefunction ipo, and 
suppose ip g is a guessed (fancy word: ansatz) wavefunction. If we define the functional 1 £ : ip H > f ip*Hip, 
then the variational theorem states that 

£[lpg\ > E 0 , 

with equality holding if and only if ip g = ipo- In other words, if ip g has optimizable parameters Ci,... ,c n , 
minimizing £[ip g (ci ,..., c n )\ with respect to these parameters will get us closer to the true ground state 
energy (and wavefunction) without having to worry about “overshooting.” 

If we have a set of normalized atomic orbitals (AOs) (pi,... , <p n for our system, we have a good place to 
start optimizing: let ip g = N(c\(p\ + ■ ■ ■ + c n <p n ), where A is a normalization constant. Thus, our molecular 
orbitals (MOs) will be formed as linear combinations of atomic orbitals (LCAOs). (Molecular orbitals 
are not really LCAOs but this assumption is a convenient approximation.) Moreover, we will assume these 
wavefunctions are all real, so that (pi = (p*. As our first order of business, we want to normalize ip g by 
determining what N should be, so we set 


J ip 2 g = N 2 J (c^ H-b c n (p n ) 2 = N 2 [c\ f (pj-\ -b c 2 n J (p 2 n + 2 ^ dCj J (pi(pj\ = 1. 

i<j 

We can simplify this by defining Sij := f (pi(pj (called the overlap integral for i yb j; Su = 1 is simply the 
normalization condition for each (pi), so we obtain the value of N: 


N = 


^ ^ C, Su -b 2 *y ' CiCjS-ij 

i i<j 



The task at hand is to minimize £ [ip g \ = J ip g Hip g : 


ipgHip g = N 2 (Ci(p! H-b c n (p n ) H [c\(p\ H-b c n (p n ) 


— N 2 [c 2 J (p\H(pi-\ - + c 2 J (p 

- .V J ;X ,: < U " + 2 V dCjHij], 

i i<j 


+ 2 ^ ^ CjCj J (j)iH(f)j\ 

i<j 

« 


where we used the abbreviation Hij := f (piH(pj (called the Coulomb integral when i = j and exchange 
integral when i yb j), and the fact (Hermitian property of the Hamiltonian) that f (piHpj = J (pjH(pi. 

When £ reaches its minimum value, the partial derivatives of £ with respect to the coefficients Cj are all 
zero: d£/dci = 0 (i = 1,... ,n). 

1 A functional is a type of function that takes a function as its input and outputs a (real) number. 
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Rewriting (*) as £ ■ [YitfSu + 2 Yi<j c i c jSij] = Yj '///,, + 2Yi<j c iCjHij and differentiating both 
sides of the equality with respect to Cj, we obtain (after combining terms) £ ■ (Yj ^ c jSij) = Yj ^ c jHij, or 
Y j Cj(H ij -£S ij ) = 0 (i = 1,2..., n). 

Using matrices, we can combine these n equations as 

Hi -| — £ S- t i H\2~£S\2 ■■■ H\ n — £S\ n ci 0 

H 21 — £S2 1 H 22 — £Sl2 • • • H2n ~ ££>2n C 2 0 

UII 1 £S n 1 H n 2 £S n 2 * * * I /1 n £Snn dl 0 

By a theorem of linear algebra, this homogeneous matrix equation has nontrivial solutions (a solution 
other than ci = • • • = c n = 0) if and only if the n x n matrix is singular: 

det {[Hij — £Sij] nxn ) = 0 (the secular determinant). 

Physical interpretation and the Hiickel assumptions 

The integrals H lt , H,j , S l3 admit relatively concrete physical interpretations. The Coulomb integral Hu 
is the energy of an electron in AO (j>i before interaction, while the exchange integral (also called the 
resonance integral) Hij (i ^ j ) is the energy of stabilization by allowing an electron to delocalize between 
<pi and (pj (or, roughly, the strength of the interaction between the orbitals). The overlap integral S t j (i ^ j) 
expresses the extent of spatial overlap between the electron densities of <t>i and <f>j . 

Hiickel MO theory is a greatly simplified theory that deals with the 7r-electrons only of a molecule and 
assumes that the cr-electrons form the molecular framework but can otherwise be ignored. In the simplest 
version of Hiickel theory dealing with unsaturated hydrocarbons, all H rl are equal and set to the energy of 
a carbon 2p orbital and is designated by a: Hu = a. 

Simplifying assumption 1: Only p orbitals on carbons i,j that are cr-bonded to each other are considered 
to interact with energy p: H l: j = /3 ( i,j neighbors). For non-neighboring and non-identical i,j, H, :l = 0. 
The signs of a and /? are negative by convention. 

Simplifying assumption 2: All overlap integrals are set to zero: Sir = 0 (i j), while because the AO’s 
are already assumed to be normalized, we can write Su = / 4 >i 4 >i = 1. 

Neglecting overlap is a drastic assumption! Moreover, there is an important consequence of having 
nonzero (positive) overlap integrals: Bonding orbitals are less stabilized than antibonding orbitals are desta¬ 
bilized: If bonding and antibonding orbitals are both filled, the overall interaction is destabilizing! However, 
we do not see this in the Hiickel approximation because overlap is set to zero. 

Example: Consider the allyl system, with atoms numbered from left to right as 1, 2, and 3. Using the 
standard Hiickel theory variable names, the secular determinant is 

a — £ p 0 

P a — £ P =0. 

0 P a — £ 

x 1 0 

To simplify notation, set x := (a — £)/p. Then equivalently, 1 x 1=0. Using minor expansion to 

0 1 x 

evaluate the determinant, we get x 3 — 2 x = 0, which has roots xi = —\/ 2 , X 2 = 0, and X 3 = \/2, which gives 
orbital energies £\ = a + V 2 p, £2 = ct, £3 — a — \/ 2 p. 

2 : 1 0 ci 0 

To work out the orbital coefficients, we plug back into the equation 1 x 1 • C 2 = 0 

0 1 x C3 0 

£y. (ci,c 2 ,c 3 ) = (c l5 v^ci.ci), so N = (c? + 2 c\ + c \)~ 1/2 = . and Vh = \<t>i + yyfo + \f>3- 

£ 2: (ci,c 2 ,c 3 ) = (ci, 0 ,-ci), so N = (cj + 0 2 + (-Ci) 2 )" 1/2 = and ^2 = ^1 - qr^- 

£3: (ci,c 2 ,c 3 ) = (d, \[2c\, Ci), so N = (cj + 2 cj + cj )~ 1/2 = and ip 3 = \<j>\ - ^</> 2 + \(t> 3 - 

Note that normalization by N = [ Yj C Y\ is equivalent to adding the condition Yi cf = 1. ■ 
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